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Abstract 

Non-compact lattice QED with two flavors of light dynamical quarks is simulated on 16 4 lattices, and 
the chiral condensate, monopole density and susceptibility and the meson masses are measured. Data 
from relatively high statistics runs at relatively small bare fermion masses of 0.005, 0.01, 0.02 and 0.03 
(lattice units) are presented. Three independent methods of data analysis indicate that the critical point 
occurs at (3 = 0.225(5) and that the monopole condensation and chiral symmetry breaking transitions 
are coincident. The monopole condensation data satisfies finite size scaling hypotheses with critical indices 
compatible with four dimensional percolation. The best chiral equation of state fit produces critical exponents 
(5 = 2.31, P m ag = 0.763) which deviate significantly from mean field expectations. Data for the ratio 
of the sigma to pion masses produces an estimate of the critical index 8 in good agreement with chiral 
condensate measurements. In the strong coupling phase the ratio of the meson masses are M^/M^ as 0.35, 
M\JM^ fa 1.4 and M%/M% « 0.0, while on the weak coupling side of the transition M%/M% « 1.0, 
M\jMp w 1.0, indicating the restoration of chiral symmetry. 
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1. Introduction and Review 



This paper continues and sharpens our series of studies of noncompact QED incorporating two species 
of light dynamical fermions. The goal of this paper will be to establish that monopole condensation in 
this theory is a legitimate phase transition with power-law critical singularities, and that this transition is 
coincident with chiral symmetry breaking. The critical properties of the chiral symmetry breaking transition 
will be studied using the theory's meson spectrum and chiral condensate. New methods of analyzing the 
simulation data of the light quark sector of the theory will suggest, as proposed in past work, that the 
critical point is characterized by a non-trivial, interacting field theory Although this paper will emphasize 
new simulation results and will leave parallel theoretical developments to other articles, we begin with some 
general remarks and perspectives. One of our motivations for this set of studies is to understand whether 
realistic strongly coupled field theories exist in four dimensions, and if so, to extract their essential ingredients. 
In the language of the renormalization groups, we are searching for ultra-violet stable fixed points which are 
non-trivial. We accept the conventional wisdom that \<p 4 is trivial in four dimensions and turn to theories 
with gauge fields. QED is the simplest candidate in this direction and if we let its charged fields be fermionic, 
chiral symmetry can protect the model from unnatural fine-tuning problems. An early excursion in this field 
found that the compact lattice version of QED coupled to dynamical fermions had a first order chiral 
symmetry breaking phase transition at strong coupling, and was not a candidate for a relativistic continuum 
quantum field theory [1] . The noncompact formulation of the model was seen to have a second order chiral 
transition which became a candidate for a continuum field theory [2] . However, since the transition occurs at 
non-zero coupling conventional classical or perturbatively-improved Actions do not necessarily yield reliable 
pictures of the theory's critical point. In fact, monopoles were discovered near the chiral transition and large 
scale simulations suggested that a topological phase transition, monopole condensation, occurs at the chiral 
transition [3]. This interesting result, which will be strengthened considerably through additional simulations 
and finite size scaling analyses here, means that conventional perturbative discussions of the triviality of QED 
are not relevant here [4] . Theories with fundamental monopoles are intrinsically nonperturbative because of 
the Dirac quantization condition. Studies of the dynamical nature of theories with fundamental electric and 
magnetic charges [5] are completely new even though the subject is fifty years old [6] . Recent simulations have 
suggested that the essential characteristics of the critical point Action driving the physics we are observing 
in simulations includes at least three ingredients: 1) a strong vector force; 2) induced four Fermi couplings; 
and 3) monopole condensation. Item 2 has been discussed and motivated elsewhere [7,8] and will not be 
discussed explicitly here. Certainly, considering the richness of items 1-3, this particular study is, alas, 
still in its infancy. Nonetheless, the article will present promising mass spectrum calculations which, when 
combined with accurate chiral condensate measurements, will set the stage for decisive critical exponent and 
anomalous dimension calculations which should help clarify the theory's critical Action. 

In Sec. 2 we introduce a random refreshment procedure to make the hybrid molecular dynamics algorithm 
more efficient. In Sec. 3 we shall present data for the chiral condensate which suggests the existence of a 
non-trivial phase transition at a coupling j3 c — 0.225(5). Since we have previously done a similar study on a 
10 4 lattice, we shall monitor possible finite size effects carefully In Sec. 4 we study monopole condensation 
in this theory and find that there is a phase transition at (3 C = 0.225 whose singularities are compatible with 
four dimensional percolation. In Sec. 5 we consider the meson mass spectrum of the model. A study of the 
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ratio of the pion and sigma masses is consistent with a power-law chiral phase transition at (3 C = 0.225(5) 
with critical indices agreeing with those found in the Equation of State study of Sec. 3. An analysis of 
the dependence of the pion mass on the chiral condensate also picks out (3 C = 0.225 and suggests that the 
transition is non-trivial. And finally, the systematics of other meson mass ratios are studied and the level 
ordering of the theory is presented. In Sec. 6 of this paper we present suggestions for further study, both 
computational and analytic. 



The lattice Action and simulation algorithm have been presented and discussed in great detail in past 
work. The reader is referred to Ref. [9] for explicit formulas. As will be discussed below, the hybrid molecular 
dynamics algorithm was considerably improved by introducing "random refreshing". This improvement 
tended to decorrelate low momentum photon modes in our gauge configurations more efficiently than in past 
studies. To control systematic errors we were forced to use a small time step in the molecular dynamics 
evolution equations. In fact, on 16 4 lattices we accumulated data with bare fermion masses m of 0.03, 0.02, 
0.01 and 0.005 in lattice units. The time step ranged from 0.02 at m = 0.03 to 0.005 at m = 0.005 to keep 
systematic errors smaller than statistical ones. The criterion for the convergence of the conjugate gradient 
algorithm used to invert the lattice Dirac operator was scaled with lattice size and fermion mass to maintain 
uniform accuracy throughout our data set. At each coupling (3 = 1/e 2 and m value we generated between 
200 and 400 molecular dynamics trajectories. Each trajectory was defined to be one time unit long, so 
400 trajectories at m = 0.005, where the time step was 0.005, corresponded to 80,000 sweeps of the hybrid 
algorithm. We collected data at couplings ranging from (3 = 0.15 (strong coupling) to (3 = 0.30 (weak 
coupling). This study dwarfs all others in the literature. Nonetheless, our meson spectrum data is not as 
accurate as we would have liked. In the last section of this paper we shall suggest further improvements on 
the present work to generate more plentiful useful results. 

In performing our updating in the weakly coupled spin wave phase we observed that the time develop- 
ment of the gauge field was very similar to that for the free gauge field, although fermion feedback shifted the 
chiral symmetry breaking transition to a considerably smaller value of the coupling (3 = 1/e 2 . The normal 
refreshment scheme is inadequate for such dynamics. We have therefore developed a random refreshment 
scheme which is adapted for the free field case and hence should be better than the standard scheme for the 
problem at hand. The free field lagrangian from which we start is 



2. A Random Refreshment Algorithm 




) 



(2.1) 



links 



sites, /i.,z/ 



In momentum space, using the Feynman gauge for convenience, L becomes 
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i.e. a set of harmonic oscillators with frequencies 
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= ^2/3(4 cos k v ) (2.3) 



so the periods Tk of these oscillators lie in the range 

-^<r fc < - 2n (2.4) 

2VP ~ ~ V 2 /?(l " cos(2tt/X)) 

For /3 = 0.225 where the transition occurs on a 16 4 lattice this range is 3.3 - 33.9. 

Ideally for any one of these oscillators, one would like to refresh often enough so that it samples its own 
(gaussian) energy distribution as rapidly as possible and loses any phase coherence with the other oscillators. 
However, we do not want to refresh it so often that it only random walks through its own configuration space. 
Our educated guess is that the best choice would be a refreshment interval ~ l/u>k [10]. If we were to choose 
this for the fastest oscillator, the slowest oscillator would random walk through its configuration space, while 
if we chose this for the slowest oscillator the fastest oscillator would spend most of its time retracing its 
orbit. We thus chose a compromise where the refreshment interval is varied randomly to cover the "optimal" 
intervals for each of the oscillators. In fact we choose the refreshment interval each time to vary between 
Tmin and r max according to 

^refresh — ^~min (j~max / ^~min} (2-5) 

where R is chosen from a set of random numbers distributed uniformly over the interval [0,1]. For our runs 
we chose r m i„ = 0.2 and r max = 10. For runs somewhat longer than ours (~ 250 time units per /3) we would 
suggest increasing T ma x closer to the maximum period of the oscillators. 



3. Chiral Equation of State 

The data for the average plaquette and < -0V > at various couplings [3 = 1/e 2 and bare fermion masses 
to are presented in Tables 1 and 2, respectively. The < -tptp > data is plotted in Fig. 1. A particularly useful, 
relatively unbiased way to extract physics from these measurements it to construct the theory's equation of 
state [11], 



< iptp > _/ A/3 



(3.1) 



m i/s y < ^ > i//3 ma3 ^ 

where A/3 = /3 — /3 C , /3 C is the bulk critical point for chiral symmetry restoration, and 5 and (3 mag are critical 
exponents. At (3 = /3 C , S records the response of the order parameter to a small symmetry breaking field, 

< ^ >= Cm 1 ' 5 (3.2) 
while Pmag controls the shape of the < > vs - P curve in the broken phase in the chiral limit m — > 0, 



< ^ >= D(/3 C - Pf^o 
4 



(3.3) 



Eq.(3.3) is difficult to use in practice because it refers to < 'tptp > at vanishing fcrmion mass m, a quantity 
which is not directly accessible to our algorithm. Our strategy will then be to exploit Eq. (3.2) to find 
candidates for (3 C and S, and then construct the equation of state Eq.(3.1) for various (3 mag . Only for the 
correct, physical values of (3 mag will all the data of Table 2 fall on the universal curve Eq. (3.1). Past studies 
on smaller lattices and polynomial extrapolations of < iptp > to zero m at fixed coupling (3, indicate that the 
critical coupling is near 0.225 and (3 mag is between 0.50 and 1.00 [3]. To search for the power-law behavior 
of Eq. (3.2) expected at the critical point, we plot — In < "tpip > vs. — lnm for m = 0.005, 0.01 and 0.02 
in Fig. 2. Note that at (3 = 0.22, the < -0-0 > value at m = 0.005 is larger than the power-law fit to the 
m = 0.01 and 0.02 data points, while at (3 = 0.225 it is smaller than the fit. This suggests that (3 = 0.220 
lies in the broken symmetry phase and the critical coupling is between 0.220 and 0.225. The straight line 
fits in Fig. 2, then predicts that the critical index S lies between 2.31 and 2.75. The reader might recall that 
this procedure was extremely successful for the quenched model where much greater statistics and smaller 
quark masses could be probed [12]. In that case data was taken with a resolution of 0.001 in f3 rather than 
0.005 used here and the critical coupling was determined to three significant digits. 

It is also interesting to search for (3 C and S in a slightly different way. Consider the plots of — 1/ ln(m) 
vs. — 1/ln < tyty > at fixed (3 in Fig. 3. This fit is particularly suitable for determining the position of the 
critical coupling and the value of the critical exponent <5. Unlike the previous plot, it brings the chiral limit, 
m = 0, to the origin. Consider the EOS in its standard form [11], m —< ipip > s f(z), with z being the 
reduced variable as defined by Eq.(3.1). The universal function f(z) is monotonically decreasing with only 
one zero in the broken phase. Denote x = — 1/ In < ^pip > and y = — 1/ ln(m). Then, 

o — x\nj(z) 

Because f(z) vanishes in the chiral limit in the broken phase, the curve y — y(x) approaches the horizontal 
axis with an infinite slope, while in the symmetric phase it approaches the origin linearly, with unit slope. 
At the critical point z = we have, 

^ = 6-l m (3 - 45) 

As we will see in later sections, the universal function f(z) will turn out to satisfy /(0) w 1 so that the 
second term in the numerator of Eq.(3.4b) can be neglected and the slope of the curve y c ~x/5 gives a fairly 
good estimate of the exponent (5. From these two equations it follows that 

--- = m^ (3.5) 
Vc y /(0) v ' 

Considering the monotonicity of the function /(z), the difference y(x) — y c {x) changes sign as the critical 
point is crossed. All the curves in the broken phase lie below, and those in the symmetric phase lie above 
the critical line. The lines lying in the broken phase terminate on the horizontal axis, whereas those in the 
symmetric phase and at the critical point terminate at the origin. Except at the critical point, all the lines 
have curvature downwards. From the figure we see that the straight line hypothesis works well at (3 = 0.225 
where it gives S — 2.29 in good agreement with Fig. 2. In addition, the curves for (3 = 0.22 and 0.210 suggest 
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that < V>V> > is non-zero in the chiral limit at m — > 0, while the curves for f3 = 0.23, 0.24 and 0.26 indicate 
that < -tpip > vanishes in the chiral limit. 

Consider the equation of state Eq.(3.1) for different hypothetical (3 C values 0.230, 0.225 and 0.220. For 
each p c and S value we can inspect Eq.(3.1) as a function of f3 mag . We will show just three characteristic plots 
from the multitude we viewed. Recall the definition of the susceptibility critical index 7, (3 mag = j/(d — 1). 
We shall show numerical evidence for the equality 7 = 1 later in this article. Recall also that 7 is exactly 
unity in the quenched model as well as in mean field theories of chiral transitions. So, we will investigate 
the hypothesis 7 = 1 and show in Figs. 4, 5 and 6, the chiral equation of state for the three cases: (1) 
(3 C = 0.220, S = 2.75, (3 mag = 0.571, (2) C = 0.225, 5 = 2.31, /3 mag = 0.763, and (3) f3 c = 0.230, 8 = 1.89, 
flmag = 1-124. The universal character of the plots - that all the data at various couplings and quark masses 
fall on one curve - is fairly good. Inspecting the three figures, we see that the dispersion in the (3 C = 0.225 
case is the least, and this fact suggests that 

Pc = 0.225, 0=2.31, (irnag = 0.763 (3.6) 

as our "preferred" conclusion from these measurements. 

These 16 4 results should be contrasted with our 10 4 results published earlier [3]. Comparing the raw 

< -0V > data between the two lattice sizes, we note considerably more precision in our present study which 
accounts for the fact that we have stronger conclusions here. Our preferred values for (3 C , S and /3 mag as 
listed in Eq.(3.6) are compatible with our 10 4 results. The only significant finite size effects in our chiral 
condensate measurements occur at the smallest m = 0.005 value where the < jpip > values on the larger 
lattice lie above those on the smaller lattice. This is the expected systematic trend - the smaller lattice 
acts as a "finite temperature" environment, suppressing chiral symmetry breaking. Note that at f} c — 0.225, 

< iptp >= 0.0881(29) on a 10 4 lattice and < ^ >= 0.0943(8) on a 16 4 lattice which is a 6% shift ignoring 
error bars and a two standard deviation effect otherwise. 

We shall see in Sec.5.b below that meson spectroscopy also favors (3 C — 0.225 and S — 2.31 in a completely 
independent fashion. 

4. Monopole Condensation 

In past work we have identified (3 C — 0.225 as the critical point for monopole condensation, using order 
parameter and susceptibilities first suggested in Ref.[13]. We can now sharpen that analysis using finite size 
scaling analyses to argue that monopole condensation is a true phase transition with power-law thermo- 
dynamic singularities. The precise definition of bond percolation and the identification of monopoles within 
a gauge field configuration can be found in Ref. [13] and will not be repeated here. Our order parameter for 
the monopole condensation transition will be M — n max /n to t which is the ratio of the number of dual sites 
in the largest monopole cluster to the total number of connected sites in all monopole clusters. Fluctuations 
in M generate a monopole susceptibility \ which should diverge at the critical point. In Fig. 7 we show \ 
vs. (3 on the 16 4 lattice for bare quark masses of 0.03, 0.02 and 0.01. There is no systematic dependence on 
the quark mass over this range of relatively small values. The critical point for monopole condensation is 
seen to be (3 C — 0.225(5). Recall that in the quenched noncompact lattice QED model x peaked at weaker 
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coupling (3 = 0.244, so charge screening has pushed the monopole condensation point to stronger coupling 
as expected [14]. However, the peak values of x are within error bars in the two cases. Finite size scaling 
relates the maximum value of x, Xmax(L), on a lattice of volume L 4 to the critical indices 7 mon and v mon of 
the monopole condensation transition, 

Xma X ~ LT-"/"""" (4.1) 

In Fig. 8 we show hiXmax plotted against InL for (3 C — 0.225 for the data presented in Table 3. For each 
lattice size we confirmed that (3 C = 0.225(5) - no size dependence was observed in the critical point of this 
transition (This peculiar fact was also observed in the quenched case [14].). Note that a linear fit to Fig. 8 is 
very good and we read off the ratio of the critical indices j mon and v moni 



Jmon/Vmon = 2.16(6) (4.2) 

In the quenched case we measured ^mon/^mon = 2.24(2). 

Now consider the behavior of the order parameter M — n max /n tot . First consider histograms of M on 
the 16 4 lattice at m = 0.01 and variable (3. In Fig. 9a we consider weak coupling (3 — 0.24 and find that M 
is essentially zero with little dispersion. However, at (3 = 0.23 shown in Fig. 9b, the distribution has moved 
to positive M and has broadened. At (3 = 0.225 shown in Fig. 9c the effect is most dramatic. Continuing to 
(3 = 0.22 (Fig.9d), (3 — 0.21 (Fig.9e) and (3 = 0.20 (Fig.9f), we see the order parameter increase in value and 
its fluctuations diminish. A plot of M vs. (3 will also be presented below. 

A particularly effective way to measure the "magnetic" critical exponent (3 mon , 

M{(3) = C(f3 c - f3f™°" , (0 < (3 C ) (4.3) 

which controls the nonanalytic vanishing of M in the L limit, is to use finite size scaling. At the critical 
point and on a L 4 lattice M((3 C , L) should vanish with a power-law, 



M(/3 C , L) - L-Pmon/Vmon ) ( 4 4 ) 

Using the data in Table 3 we show In M((3 C ,L) vs. InL in Fig. 10 for L ranging from 10 to 18. The error 
bars are much larger than we would prefer, but an estimate of (3 mon /v mon) follows from the L = 12, 14, 16 
and 18 points: 



Pmon/Vmon = 0.64(22) (4.5) 

Similar measurements in the quenched model [14] gave (3 mon /v mon = 0.88(2). The values of M((3 C ,L) for 
the quenched model were systematically below those of Table 3. It is not clear whether fcrmion feedback 
has reduced l3 mon /v mon , below its quenched value since the error bars in Fig. 10 are relatively large. 

In Ref.[14] we argued that monopole condensation in quenched noncompact lattice QED is a percolation 
process uncomplicated by dynamical considerations. Our perspective on this problem comes from bond 
percolation models where one randomly occupies bonds with a probability p. If p is small, typical lattice 
configurations have small, improbable connected clusters of occupied bonds. However, as p increases towards 
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a critical concentration p c , the largest connected cluster becomes infinite in extent and occupies a macroscopic 
fraction of the available lattice bonds. The value p c is called the percolation threshold and traditional 
statistical mechanics calculations produce the estimate p c w 0.16 in four dimensions. Following Rcf.[13] we 
computed the total amount of monopole current £, £ — |m^(n)|, where m^(n) is the integer-valued 

monopole current defined on links of the dual lattice, in each of our 16 4 configurations. The resulting 
curve of the average value of £ normalized by the number of lattice links is shown in Fig. 1 1 . We also 
show the percolation order parameter M = n max /n tot there for the same set of configurations having small 
m = 0.01. We see that M turns on at (3 — 0.225 with the expected finite size rounding and that < £ > /Ng 
is roughly 0.135 and varies smoothly around the critical point. We learn that < I > /N t is near the critical 
concentration of pure bond percolation and since it varies smoothly it could be used as an independent 
variable to parameterize the monopole measurements instead of the gauge couplings (3. It is interesting that 
< I > /Ni varies abruptly across the confinement transition of pure compact U(l) lattice gauge theory [15]. 
The character of Fig. 11 suggests that the monopole dynamics of the two flavor noncompact model is not 
qualitatively different from the quenched model where accurate finite size scaling measurements placed the 
monopole transition in the universality class of four dimensional percolation [14]. 

In summary, the two flavor theory has a monopole condensation transition at (3 C — 0.225 with power- 
law thermodynamic singularities. The determination of 7 mo „/^ mo „ is in good agreement with the same 
calculation for the quenched model which in turn has critical indices in agreement with four dimensional 
percolation [16]. The measurement of (3 m on/^mon is less quantitative, but the data for M — n max jn tot cer- 
tainly satisfies a scaling hypothesis. It is tempting to conclude that the chiral symmetry breaking transition 
and the monopole condensation transition are coincident. If that is true, then the set of field theory critical 
indices controlling the chiral transition (<5 = 2.31, (3 mag = 0.763) would be distinctly different from mean 
field values. 

It is sometimes said that the monopoles of lattice QED are artifacts of the regularization scheme and 
have no physical significance. However, in this case we have seen that they experience a second order phase 
transition with scaling laws and critical indices which satisfy hyperscaling. Therefore, their critical properties 
have decoupled from the lattice itself and they should have macroscopic physical consequences. 



5. Spectrum 

a. Raw Data 

We calculated the propagators of the pion, sigma, rho, and A\ mesons using standard techniques bor- 
rowed from lattice QCD large scale simulations [17]. Most of the meson propagator calculations used a point 
source for the quark propagator calculation and then took the zero momentum projection of the propagator 
perpendicular to the "temporal" direction of the 16 4 lattice. The fall-off of the meson propagator in the 
temporal direction then determined the mass of the state using fitting technology borrowed from lattice 
QCD [17]. The standard parameterization reads, 

G(t) = Q (exp(-Mir) + exp(-M i (16 - r))^ + ^](-l) T A ^cx P (-M 4 t) + exp(-M,(16 - r))^ (5.1) 
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and we considered fits over different ranges of r between and 8. fn most cases good fits were found with just 
the leading single particle term in Eq.(5.1) determining the amplitude C, and the mass M, but in several 
cases secondary particles were also necessary although their masses and amplitudes were not determined 
with sufficient accuracy to be interesting. Error estimates on the masses followed from the usual correlated 
X 2 procedures, also borrowed from lattice QCD studies [17]. 

Before turning to the masses themselves, it is instructive to consider typical fits. Consider the coupling 
(3 = 0.225 and m = 0.01, and the pion (Fig. 12), sigma (Fig. 13), rho (Fig. 14) and A\ (Fig. 15) propagators. 
We show both the raw data (open circles) and the fits (triangles). As usual, the fits for the pion were the 
best, and we will be able to do some quantitative analysis with the data here. The sigma particle also 
admitted respectable fits and its mass estimates will be useful. The rho and A\ were relatively disappointing 
and very little quantitative information will be obtained in these cases. 

The data of meson masses (squared) for couplings ranging from 0.20 to 0.26 and quark masses from 
0.03 to 0.005 are given in Tables 4, 5, 6 and 7. The quality of the fits is indicated in part by the error bars 
accompanying each entry. Typically, especially in the broken symmetry phase, the confidence levels of the 
fits were very good (90% and higher confidence levels were common). However, the quality of the fits in the 
weak coupling phase was distinctly poorer. The simple hypothesis Eq.(5.1) is probably too naive in the weak 
coupling phase where light charged fermions propagate freely, interacting with massless photons, etc. This 
result is certainly not very surprising, but we did not know better fitting forms than Eq.(5.1) which could 
be used productively on our limited data set. In addition, examining the rho and A\ tables, we see that the 
large statistical errors in the data renders the fits of little use. Missing entries in the tables indicate that the 
fitting routine was unable to find an acceptable fit for that particular propagator. Clearly it is meaningless 
to extrapolate the rho and A\ masses to m — > 0. Therefore, our analysis will consider m = 0.01 for the most 
part for these two states. Happily, the pion and sigma results are much more quantitative and we will learn 
a lot from their m and f3 dependencies. 

b. The Goldstone Pion 

If chiral symmetry is really broken at (3 C — 0.225, then the pion should be massless in the strong 
coupling phase. This effect is quite clear in the data. If we extrapolate the pion mass-squared to m = 0, 
zero symmetry breaking, then Fig. 16 follows. The m — > extrapolation was done using the m = 0.03, 0.02 
and 0.01 data and the form 

Ml = a + bm + cm 2 + ... (5.2) 

expected from PCAC. Although the error bars in Fig. 16 are quite large, M% is compatible with zero for 
{3 < 0.225 and is non-zero in the weak coupling phase. Recall, however, our caveat about the fitting procedure 
Eq.(5.1) at weak coupling the simple pole structure of composite meson propagators in momentum space 
is probably too naive in the weak coupling phase to parameterize the physics quantitatively. Nonetheless, 
the qualitative fact that the pion mass is distinctly non-zero in the weakly coupled phase is confirmed by 
the simulation. Chiral symmetry restoration in the weak coupling phase will be discussed in greater detail 
below. 
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c. Spectroscopy and the Critical Index 6 

We have argued in Ref.[18] that the properties of ratios of dynamically generated masses M a and Mb, 
R a b(m, 0) = G(m/ {0 — /3 C ) A )), with A = S0 ma g, are determined by the Equation of State of the chiral order 
parameter. Let R(m, (3) denote the ratio of M% to M 2 and write the Equation of State in the standard form, 

A/3 

< > 1/A 



<^>°/( ; T A t a ) (5.3) 



with A0 = C - (3. Then [18], 



1 -x * /'(») (5.4a) 



R(m,f3) I3 mag f{x) 

where x is the reduced coupling, A/3/ < -tpip > 1 // 3 ™»s. Eq.(5.4a) reduces at the critical point — C (x = 0) 
to, 

M 2 1 

This result will give us another measurement of C and 6 to supplement and, in fact, support Eq. (3.6). 
Recall that the derivation of Eq. (5.4a-b) follows from the chiral equation of state describing the spontaneous 
symmetry breaking of a continuous group [18]. In fact, the derivation reveals more than Eq. (5.4b) which 
states that at the critical point the ratio of the pion and sigma masses squared is independent of m and is 
equal to 1/5. For fixed < C in the broken symmetry phase, i?(m, 0) should fall to zero as m — > since 
the pion is a Goldstonc boson. For fixed (3 > (3 C in the symmetric phase, R{m, (3) should increase to unity 
as m — > since the sigma and pion must be degenerate in the symmetric phase. Thus, in the scaling region, 
the following bound can be placed on S [18] t 

R(m,0<0 c )<±<R(m,0>0 c ). (5.5) 
o 

Because of the scaling form of the mass ratio, this bound improves as either approaches C or as m becomes 
large. 

In Fig.17 we plot R(m,f3) using the data for M% and M 2 in Tables 4 and 5. For (3 = 0.20, 0.21 and 
0.22 R clearly falls as m decreases while for (3 = 0.23, 0.24 and 0.25 R increases as m decreases. Only at 
(3 = 0.225 is R independent of m and has the value 0.391(11) which gives 



S = 2.56(7) (5.6) 

The consistency of these results with the analysis based on the chiral equation of state in Sec. 3 is quite 
striking. In fact, the determination of C = 0.225(5) is most persuasive from Fig.17 - R(m, 0) systematically 
decreases with m at = 0.22 and systematically increases as m decreases at = 0.23, thus picking out 
0c — 0.225(5) as the critical point. 

t An analysis very similar to ours was used to determine the universality class of the finite temperature 
chiral restoration transition in QCD at infinite couplings by, G. Boyd, J Fingberg, F. Karsch, L Karkkainen 
and B. Peterson, Nucl. Phys. B376, 199 (1992). 
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It is interesting to pursue Eq. (5.3) and (5.4) in somewhat greater detail. In particular, Eq. (5.4a) 
states that the mass ratio R = M%/M% follows from the Equation of State with no additional parameters, 
etc. Do our chiral condensate measurements produce a universal Equation of State f(x) which predicts the 
correct systematics for R? This is a highly nontrivial test of the hypothesis that there is a chiral critical point 
at C = 0.225 with S = 2.31 and /3 ma g = 0.763 since the spectroscopy measurements of Ml and M 2 were 
independent of chiral condensate considerations. In Fig. 18 we plot the Equation of State in the "standard" 
form Eq. (5.3). One of the advantages of plotting the Equation of State in this form is that the universal 
function / is so simple. Clearly a linear fit for / represents the data adequately. In fact, as discussed in 
Rcf.[18] a linear / occurs only in a theory where the susceptibility index 7 is unity. As mentioned in our 
discussion of Fig. 4-6 and now seen in Fig. 18, the simulation data does provide direct evidence for 7 = 1. 
The interested reader is directed to Ref.[18] for further detail, since here our emphasis is on spectroscopy 
and the exponent S. A linear fit to f(x) gives, 

f(x) = -5.3125a; +1.075 (5.7) 

Using Eq.(5.4a) R(m, 0) is now determined in a universal form and is plotted in Fig. 19. Comparing Fig. 17 and 
Fig. 19 we find good agreement within the rather substantial error bars for our spectroscopy measurements. 
Nonetheless, the general success of Eq.(5.3) and (5.4) to describe spectroscopy data within the hypothesis 
of a chiral critical point with power-law non-mean field singularities is encouraging. 

d. M^ vs. < -0-0 > 2 , the Critical Point and Anomalous Dimensions 

As we have discussed in Ref.[18], measurements of the pion mass and the chiral order parameter can 
be exploited in unison to locate critical points and measure anomalous dimensions. The two theoretical 
ingredients underlying this strategy is 1) correlation scaling; and 2) the equation of state. In addition, the 
Goldstone nature of the pion, which we have confirmed explicitly above, allows additional information to be 
extracted from plots of Ml vs. < ipip > 2 - We refer the interested reader to Ref.[18] for derivations. Here 
it is more appropriate to quote results and apply them to our simulation data. It is convenient to express 
the results in terms of r], the anomalous dimension of < tyty >. According to hyperscaling the anomalous 
dimension i] is related to the critical indices (3 mag of the order parameter and v of the correlation length by 
fimag/v = 1 + rj/2. Ref.[18] predicts the squared pion mass' dependence on the square of the chiral order 
parameter within the scaling region for fixed coupling 0. Consider the broken symmetry phase < C , the 
critical point — C and the symmetric phase < C in turn, 

< C : Ml ~< > 2 - < ^ > (5.8a) 
= C : Ml~ (< <H> > 2 ) 1 /(i+V2) ( 5 86 ) 

> C : Ml~(<W >2)l/(l+2/r,) (5 8c) 

So, Eq.(5.8a) states that Ml should vanish linearly asm^O and < tpip > 2 ^>< iptp >g^ 0. Eq.(5.8b) states 
that at the critical point Ml vs < -0V > 2 should be concave downward and provide a measurement of the 
anomalous dimension r\. Recall that a physically meaningful theory must have r\ > 0, so 1/(1 + rj/2) < 1. 
Finally, Eq.(5.8c) states that the curves in the symmetric phase will also be concave downward and will be 
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sensitive to 77. Clearly, Eq.(5.8c) docs not apply in the case 77 = 0. The more general formulas of Ref.[18] 
show that in a mean field scenario where 77 = that M 2 should vary as C(0—0 C )+ < > 2 m the symmetric 
phase. 

In summary, we learn from Eq.(5.8) that plots of M 2 vs. < ipip > 2 are sensitive probes into the 
dynamics of the critical point. For example, if mean field theory applies then the M 2 curves should be linear 
functions of < ■ipip > 2 for all coupling and the curve which passes through the origin picks out the critical 
point uniquely. If an interacting field theory describes the transition, then the M 2 curves should be linear 
only in the broken symmetry phase and otherwise have downward curvature determined by the anomalous 
dimension 77. 

In Fig. 20 we present our M 2 and < -tpip > 2 data in this format. Clearly the curves are linear for 
< 0.220 and give a non-zero chiral condensate in the chiral limit M 2 — > 0. The functional dependence of 
the = 0.225 and 0.230 curves is not decisively predicted, while the = 0.240 and 0.250 curves are concave 
downward. We obtain an estimate of C = 0.225 — 0.230 in good agreement with our other approaches and, 
although our data is not good enough to yield a useful value for 77, wc obtain qualitative evidence that the 
critical point is not described by a mean field theory. 

e. Mass Ratios and Chiral Symmetry Restoration at Weak Coupling 

One of our aims in these simulations is the determination of meson mass ratios in the chiral limit. As 
shown in our discussion of M 2 /M 2 these quantities can indicate whether the chiral critical point in trivial 
and teach us about the realization of chiral symmetry in the two phases. Unfortunately our rho and A\ data 
are not quantitative enough to use in a similar fashion. Therefore, we shall only consider m = 0.01 data 
here and look for trends which we believe should survive the m — > extrapolation. 

First consider the ratios M 2 /M 2 and M 2 /M 2 at variable coupling 0. In Fig. 21 we present the results 
which indicate that both ratios have non-zero values in the weak coupling phase. In fact M 2 /M 2 rapidly 
increases to 0.915(15) at = 0.26. We expect M 2 /M 2 = 1 in a phase where chiral symmetry is not 
spontaneously broken, so this is a welcome result. Looking back at Fig. 17 we see R(m, 0) rapidly rising 
toward unity as m — > for values of 0.24 and 0.25, so this particular result looks rather firm. Note from 
Fig. 21 that M 2 /M 2 appears to be non-zero in the weak coupling phase but its particular value is not related 
to symmetry considerations. 

Next, consider the ratio Mp/M\ i at m = 0.01 for various couplings as plotted in Fig. 22. At weak 
coupling, = 0.26 say, the ratio is consistent with unity, which is the expected value in a phase where chiral 
symmetry is an unbroken symmetry, 

= 0.96(8) (weak coupling) (5.9) 

As we pass to the strong coupling phase this ratio has certainly decreased, but its precise value is obscured 
by large error bars. The data suggests considerable splitting in the strong coupling phase, 

Ml 

w 1.35 (broken phase) (5.10) 
but we hesitate to place an error bar on Eq.(5.10). 
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The ratio M% /Mp 1 at m — 0.01 varies much less dramatically as /3 is varied. In the vicinity of the critical 
point and at weak coupling, Fig. 23 shows that 




/M^ is plotted against 



0.35 



/M\ 



r p but the figure is obscured 



(5.11) 



6. Conclusions and Prospects 



Instead of summarizing our results, we shall end with a few words about future work. It should be clear 
to the reader, however, that considerable progress has been made in understanding the two flavor version of 
noncompact lattice QED and its chiral transition. That transition appears to be characterized by non-mean 
field exponents and monopoles probably play an important role in its existence. The use of the Equation 
of State and spectroscopy data appears to be a particularly productive approach towards establishing its 
physical characteristics. 

Obviously we must confirm these results on larger lattices and push our finite size scaling studies further 
before solid claims can be made about the continuum limit of the lattice theory. Spectroscopy studies using 
wall sources should be made on asymmetric lattices with ten times the statistics accumulated here. A serious 
simulation on 64 x 32 3 lattices seems essential. 

The physical basis of the chiral transition at (3 C = 0.225(5) must be clarified. We need to obtain a quan- 
titative understanding of the importance of monopoles at the transition. We have already begun a project 
in that direction. By cooling lattice configurations already generated here, we can produce configurations 
with the monopoles intact, but the photons removed. Is < -0-0 > the same in these cooled configurations 
for (3 < 0.225? Additional calculations of this variety may clarify the content of the lattice theory's critical 
Action. Interesting work in this area has already appeared [19]. What is the fate of the fermion and the 
photon in the broken symmetry phase? Our attempts to do direct measurements of these propagators were 
not informative due to statistical noise. In the future we hope to address these crucial questions using 
improved wavefunction techniques borrowed from lattice QCD spectrum calculations. 



The computer simulation data presented and analyzed here required 3000 hours on one 8K quad of the 
Connection Machine (CM-2) of NCSA. Simulations were also done using the CM-2 at NPAC and PSC, while 
the ST-100 of Argonne National Laboratory contributed some data before its untimely demise. The work of 
J. B. Kogut is supported in part by the National Science Foundation, NSF-PHY97-00148. D. K. Sinclair is 
supported by DOE contract W-31-109-ENG-38. 
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Table 1 



Average plaquette data at various bare fermion masses to and couplings (3 — 1/e 2 on a 16 4 lattice. 






to = .005 


TO = .010 


to = .020 


to = .030 


.30 






0.7803(2) 


0.7815(2) 


.28 






0.8332(2) 


0.8338(2) 


.27 






0.8621(3) 


0.8625(3) 


.26 




0.8922(3) 


0.8943(2) 


0.8949(3) 


.25 




0.9264(4) 


0.9283(2) 


0.9296(3) 


.24 




0.9631(4) 


0.9651(3) 


0.9662(3) 


235 






0.9858(3) 


0.9869(3) 


.23 


1.0025(2) 


1.0043(4) 


1.0070(3) 


1.0076(2) 


225 


1.0230(4) 


1.0266(4) 


1.0290(3) 


1.0302(3) 


.22 


1.0465(3) 


1.0489(4) 


1.0521(3) 


1.0539(3) 








l.Ui (£y±) 


l.U <oo(o j 


.21 




1.1011(4) 


1.1049(5) 


1.1047(3) 


.20 




1.1607(9) 


1.1635(5) 


1.1628(4) 


.19 






1.2300(5) 


1.2276(4) 


.18 






1.3028(7) 


1.2988(4) 


.17 






1.3883(6) 


1.3818(5) 


.16 






1.4821(7) 


1.4757(4) 


.15 






1.5873(9) 


1.5812(4) 
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Table 2 



< tpip > data at various bare fermion masses to and couplings /3 = 1/e 2 on a 16 4 lattice. 



(3 to = .005 m = .010 m = .020 to = .030 



.30 






.0567(1) 


.0821(1) 


.28 






.0697(1) 


.0988(1) 


.27 






.0793(2) 


.1111(2) 


.26 




.0499(2) 


.0924(2) 


.1265(2) 


.25 




.0630(4) 


.1098(3) 


.1452(3) 


.24 




.0824(5) 


.1331(3) 


.1675(3) 


.235 






.1483(4) 


.1812(2) 


.23 


.0764(9) 


.1135(8) 


.1636(6) 


.1955(2) 


.225 


.0943(8) 


.1340(7) 


.1809(5) 


.2119(3) 


.22 


.1225(9) 


.1539(6) 


.1980(6) 


.2288(4) 


.215 






.2200(11) 


.2462(2) 


.21 




.2056(9) 


.2395(9) 


.2638(3) 


.20 




.2585(14) 


.2827(8) 


.3019(4) 


.19 






.3261(7) 


.3409(4) 


.18 






.3678(11) 


.3772(3) 


.17 






.4085(10) 


.4138(3) 


.16 






.4427(7) 


.4481(2) 


.15 






.4774(9) 


.4792(2) 
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Table 3 



Peak of the monopole susceptibility and the order parameter M = n max /n to t at (3 C = .225 on L 4 lattices 
with L ranging from 10 to 18. 



L 


Xmax 


M 


10 


53.3(2.4) 


.222(12) 


12 


78.3(4.1) 


.212(13) 


14 


108.7(6.2) 


.193(14) 


16 


143.1(5.2) 


.177(7) 


18 


190.4(8.8) 


.1630(53) 
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Table 4 



values for various couplings (3 and bare fermion masses m. 



m\(3 


.20 


.21 


.22 


.225 


.23 


.24 


.25 


.26 


.03 


0.229(3) 


0.246(3) 


0.266(3) 


0.264(12) 


0.279(14) 


0.272(45) 


0.287(27) 


0.293(25) 


.02 


0.159(2) 


0.176(2) 


0.197(2) 


0.212(4) 


0.207(4) 


0.237(5) 


0.254(5) 


0.296(6) 


.01 


0.083(2) 


0.096(3) 


0.114(4) 


0.124(3) 


0.137(3) 


0.171(9) 


0.214(5) 


0.250(5) 


.005 






0.069(3) 


0.082(2) 


0.095(2) 









Table 5 

values for various couplings (i and bare fermion masses to. 



m\/3 .20 .21 .22 .225 .23 .24 .25 .26 

.03 1.051(49) 0.992(32) 0.709(22) 0.681(16) 0.662(11) 0.567(17) 0.505(10) 0.522(8) 

.02 1.136(78) 0.716(48) 0.570(16) 0.520(20) 0.455(13) 0.452(15) 0.380(8) 0.386(7) 

.01 0.610(87) 0.475(22) 0.358(12) 0.339(11) 0.286(9) 0.246(7) 0.252(5) 0.273(5) 

.005 0.233(11) 0.205(6) 0.173(6) 



18 



Table 6 

m 2 p values for various couplings f3 and bare fermion masses to. 



m\(3 .20 .21 .22 .225 .23 .24 .25 .26 

.03 2.17(15) 1.073(93) 0.412(54) 1.36(11) 1.24(11) 1.195(71) 0.789(34) 0.837(35) 

.02 0.83(32) 1.25(91) 1.17(40) 0.645(52) 0.564(99) 

.01 1.15(25) 1.11(18) 1.03(17) 0.956(88) 0.817(75) 0.776(62) 0.767(55) 0.663(33) 

.005 0.56(19) 0.98(32) 0.57(11) 



Table 7 

m 2 A values for various couplings and bare fermion masses to. 



m\(3 .20 .21 

.03 2.19(1.23) 

.02 1.80(1.23) 1.99(80) 

.01 0.87(48) 0.74(28) 
.005 



.22 .225 .23 

1.19(23) 1.69(38) 1.61(29) 

1.46(28) 1.70(21) 1.85(19) 

1.39(35) 1.08(15) 1.18(14) 



.24 .25 .26 

1.49(23) 1.113(97) 1.201(85) 

1.43(15) 1.580(97) 1.160(72) 

0.920(94) 0.808(62) 0.693(35) 
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Figure Captions 

1. < ifnp > vs. (3 for m = 0.03, 0.02, 0.01 and 0.005 for 16 4 lattice. The data is given in Table 1. The 
size of the symbols in the figure includes the statistical error bars. 

2. - In < -0^ > vs. - In m plots for /3 = 0.22, 0.225 and 0.23. 

3. — l/ln(m) vs. — 1/ln < tpip > plots for various (3 values. 

4. Equation of State Eq.(2.1) for (3 C = 0.220, 6 = 2.75, (3 mag = 0.571. 

5. Equation of State Eq.(2.1) for /3 C = 0.225, S = 2.31, /3 mag = 0.763. 

6. Equation of State Eq.(2.1) for (3 C = 0.230, S = 1.89, f3 mag = 1.124. 

7. Monopole susceptibility \ vs - /? on a 16 4 lattice and bare fermion masses 0.03, 0.02 and 0.01. 

8. In 

Xmax vs - l n L plotted from the data in Table 2. The lattices range in linear dimension L — 10, 12, 14, 
16, and 18, and the coupling is (3 C = 0.225. 

9. a. Histogram of the monopole order parameter at (3 = 0.24. 

b. Histogram of the monopole order parameter at (3 = 0.23. 

c. Histogram of the monopole order parameter at (3 = 0.225. 

d. Histogram of the monopole order parameter at (3 = 0.22. 
c. Histogram of the monopole order parameter at [3 — 0.21. 
f. Histogram of the monopole order parameter at (3 = 0.20. 

10. lnM(/3 c , L) vs. ln L from the data in Table 3 at (3 C = 0.225. 

11. The monopole concentration < t > /Ni (triangles) and the percolation order parameter M (circles) on 
a 16 4 lattice at m = 0.01. 

12. A typical pion propagator and its fit, (3 = 0.225 and m = 0.01. 

13. A typical sigma propagator and its fit, (3 = 0.225 and m = 0.01. 

14. A typical rho propagator and its fit, (3 = 0.225 and m = 0.01. 

15. A typical A\ propagator and its fit, (3 = 0.225 and m = 0.01. 

16. M% vs. (3 at to — > 0. The plot was obtained by extrapolating the data in Table 4 to to = as discussed 
in the text. 
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17. R(m,0) vs. m for various couplings 0.200 (circles), 0.210 (triangles), 0.220 (inverted triangles), 0.225 
(hexagons), 0.230 (dark circles), 0.240 (dark squares) and 0.250 (dark triangles). 

18. The chiral Equation of State plotted in the form Eq.(5.3) with C = 0.225, 6 = 2.31 and (3 mag = 0.763. 
The line is the best linear fit to the universal function /. 

19. R(m,f3) vs. m from Eq.(5.4a) and (5.6) at values coinciding with Fig. 17. 

20. Ml vs. < Tpip > 2 for the data of Tables 2 and 4. 

21. Ml/Ml and M 2 /M 2 vs. for fixed m = 0.01. 

22. M 2 p /M 2 Ai vs. (3 for m = 0.01. 

23. M ff 2 /M p 2 vs. for m = 0.01. 

24. Mj/M 2 vs. M 2 /M p 2 for m = 0.01. 
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